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In the theory of quantum dynamical filtering, one of the biggest issues is that the underlying 
system dynamics represented by a quantum stochastic differential equation must be known exactly 
in order that the corresponding filter provides an optimal performance; however, this assumption 
is generally unrealistic. Therefore, in this paper, we consider a class of linear quantum systems 
subjected to time-varying norm-bounded parametric uncertainties and then propose a robust ob- 
server such that the variance of the estimation error is guaranteed to be within a certain bound. 
Although in the linear case much of classical control theory can be applied to quantum systems, 
the quantum robust observer obtained in this paper does not have a classical analogue due to the 
system's specific structure with respect to the uncertainties. Moreover, by considering a typical 
quantum control problem, we show that the proposed robust observer is fairly robust against a 
parametric uncertainty of the system even when the other estimators — the optimal Kalman filter 
and risk-sensitive observer — fail in the estimation. 
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I. INTRODUCTION 

Quantum filtering theory was pioneered by Belavkin 
in remarkable papers [l|, 0, 9 and was more lucidly re- 
considered by Bouten et al. 0, [1[ . This theory is now 
recognized as a very important basis for the development 
of various engineering applications of the quantum the- 
ory such as quantum feedback control [f| 0, 1MJ1, E 
quantum dynamical parameter estimati on J 12 . \\l 
and quantum information processing [lH Il6j |. 

We here provide a brief summary of the quantum fil- 
tering theory by using the same notations as those in 
0, [E[ . Let us consider an open system in contact with a 
field, particularly a vacuum electromagnetic field. This 
interaction is completely described by a unitary operator 
tit that obeys the following quantum stochastic differen- 
tial equation (QSDE) termed the Hudson-Parthasarathy 
equation [17] : 



hdU t = 



1 



-iff - -c r c)dt + cdBj 



c)dB t 



U u U = I, 
(1) 

where c and ff are the system operator and Hamiltonian, 
respectively. The quantum Wiener process B t , which is 
a field operator, satisfies the following quantum Ito rule: 

dB t dB t = 0, dB[dB t = 0, dB t dB f t = hdt, dB\dB\ = 0. 

The time evolution of any system observable X under 
the interaction (fTJ is described by the unitary transfor- 
mation jt(X) := \j\XUt- The infinitesimal change in 
this transformation is calculated as 

Mj t (X) =j t (LX)dt+j t ([tf,X})dB t +j t ([X,c})dBl (2) 



Here, we have defined CX := i[ff , X] + c) Xc - \c)cX - 
hX&c. The field operator after the interaction is de- 
termined by B' t :— j t (B t ). In the homodyne detec- 
tion scheme, we measure the field operator of the form 
Y t := B' t + B'J, which results in 



dY t = j t (c + tf)dt + dB t + dB 



(3) 
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An important fact is that the above observable is self- 
nondemolition: [Y Sl Y t ] — for all s and t. This im- 
plies that the observation is a classical stochastic pro- 
cess. (For this reason, we omit the "hat" on Y t , but 
note that it itself is not a c-number.) It is also notewor- 
thy that Y t satisfies the quantum nondemolition ( QND ) 
condition, [Y s ,j t (X)] — Vs < t, for all system observ- 
ablcs X. Our goal is to obtain the best estimate of 
the system observable jt(X) based on the observations 
Y s (0 < s < t), which generate the Von Neumann al- 
gebra yt — vN(F s : < s < t). As in the case of the 
classical filtering theory, the best estimate in the sense of 
the least mean square error, ((j t (X) — X') 2 ) — > min., is 
given by (a version of) the quantum conditional expec- 
tation: X' = Tit (X) := V(j t (X)\y t ). Here, the expec- 
tation (X) is defined by (X) := Tr [X(p ® $)], where p 
and <& represent the system quantum state and the field 
vacuum state, respectively. It should be noted that the 
following two conditions must hold in order for the above 
quantum conditional expectation to be defined: First, 
is a commutative algebra, and second, jt(X) is included 
in the commutant of y t . But these conditions are actu- 
ally satisfied as shown above. Consequently, the optimal 
filter for the system dynamics is given by the change 
in nt(X) as follows: 

hdTTt(X) = n t (CX)dt + [ir t {Xc+ c f A) 

- 7T t (A> t (c + gt)] [ d Y t - 7r t (c + gt)tft] . (4) 

We can further incorporate some control terms into the 
above equation. Typically, a bounded real scalar control 
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input u t , which should be a function of the observations 
Y s up to time t, is included in the coefficients of the 
Hamiltonian. We lastly remark that the conditional sys- 
tem state pt is associated with the system observable by 
the relation nt(X) — Tr (Xpt), which leads to the dy- 
namics of pt termed the stochastic master equation. 

A key assumption in the filtering theory is that per- 
fect knowledge about the system dynamics model ([2]) is 
required in order that the filter provides the best es- 
timate of the (controlled) system observable. However, 
this assumption is generally unrealistic, and we depend 
on only an approximate model of the system. This not 
only violates the optimality of the estimation but also 
possibly leads to the instability of the estimation error 
dynamics. This problem is well recognized in the clas- 
sical filtering theory and various alternative estimators 
for uncertain systems, which are not necessarily optimal 
but robust to the uncertainty, have been proposed. (We 
use the term "filter" to refer to only the optimal esti- 
mator.) For example, in a risk-sensitive control prob- 
lem in which an exponential-of-integral cost function is 
minimized with respect to the control input, it is known 
that the corresponding risk-sensitive observer enjoys en- 
hanced robustness property to a certain type of system 
uncertainty [H, [H, 0|. Moreover, by focusing on spe- 
cific uncertain systems, it is possible to design a robust 
observer such that the variance of the estimation error is 
guaranteed to be within a certain bound for all admissi- 
ble uncertainties [Hi [H, HI, HI • 

It is considered that the above mentioned robust esti- 
mation methods are very useful in the quantum case since 
it is difficult to specify the exact parameters of a quantum 
system in any realistic situation: for instance, the total 
spin number of a spin ensemble With this back- 

ground, James has developed a quantum version of the 
risk-sensitive observer for both continuous (2f| and dis- 
crete cases [26| and applied it to design an optimal risk- 
sensitive controller for a single-spin system. We should 
remark that, however, the above papers did not provide 
an example of a physical system such that the quantum 
risk-sensitive observer is actually more robust than the 
nominal optimal filter. 

Therefore, in this paper, we focus on the robust ob- 
server and develop its quantum version. More specifi- 
cally, we consider a class of quantum linear systems sub- 
jected to time-varying norm-bounded parametric uncer- 
tainties and obtain a quantum robust observer that guar- 
antees a fixed upper bound on the variance of the estima- 
tion error. Although in the linear case much of classical 
control theory applies to quantum systems, the robust 
observer obtained in this paper does not have a classical 
analogue in the following sense. First, unlike the classical 
case, the error covariance matrix must be symmetrized 
because of the noncommutativity of the measured system 
observables. Second, due to the unitarity of quantum 
evolution, the uncertainties are included in the system 
representation in a different and more complicated way 
than those in the classical system considered previously; 



as a result, both the structure of the quantum robust ob- 
server and the proof to derive it differ substantially from 
those found in [2l|, [H, [H Hi] . The other contribution of 
this paper is that it actually provides a quantum system 
such that both the robust observer and the risk-sensitive 
observer show better performance in the estimation error 
than the nominal optimal filter. 

This paper is organized as follows. Section II provides 
a basic description of general linear quantum systems, 
in which case the optimal filter Q is termed the quan- 
tum Kalman filter. In addition, we derive a linear risk- 
sensitive observer. In both cases, an explicit form of the 
optimal control input is provided. The quantum version 
of the robust observer is provided in Section III. Section 
IV discusses robustness properties of the proposed robust 
observer and the risk-sensitive observer by considering 
a typical quantum control problem — feedback cooling of 
particle motion. Section V concludes the paper. 

We use the following notations: for a matrix A = (ay ) , 
the symbols A T and A* represent its transpose and ele- 
mentwise complex conjugate of A, i.e., A T = (a^) and 
A* = (a*j) = (A^) T , respectively; these rules can be ap- 
plied to any rectangular matrix including column and row 
vectors. A Hermitian matrix A — A^ is positive semidef- 
inite if Av > for any vector v; the inequality A > B 
represents the positive semidefiniteness of A — B. 

II. LINEAR QUANTUM SYSTEM 

A. Quantum Kalman filter 

In this paper, we consider a single one-dimensional 
particle interacting with a vacuum electromagnetic field. 
The extension to the multi-particle case is straightfor- 
ward 11]. In particular, we focus on the particle posi- 
tion q and momentum p. The system Hamiltonian and 
operator are respectively given by 

H = ^x T Gx - x T ZBu t , c = Cx, (5) 

where x — [q p] T . Here, u t € K is the control input, 
B 6 M 2 is a column vector, and C € C 2 is a row vector. 
The 2x2 matrix G is real symmetric and £ is given by 



Then, by defining x t = [q t Pt] T = [jt{q) jt(p)] T and not- 
ing the commutation relation [q, p] = ih, the system 
dynamics ((2]) leads to the following linear QSDE: 

dx t = Ax t dt + Bu t dt + iE[C T dBl - &dB t ], (6) 

where the matrix A is defined by A :— £[G + lm(C'C)]. 
The output equation ([3]) becomes 

dY t = Fxtdt + dBt +dB\, F:=C + C*. 
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It follows from Eq. (|3]) that the best estimate of the 
system observable, n t (x) :— [7it(g) 7r f (p)] T S K 2 , obeys 
the following filter equation: 

di7 t (x) — An t (x)dt + Bu t dt 

^V t F T + £ T Im(C*) T ] (dY t - Fn t {x)dt). (7) 

In the above equation, Vt represents the symmetrized 
covariance matrix defined by 



Vt 
P, 



v{p t \y t ) 



A<?1 



UAq t Ap t + Ap t Aq t 



\{Aq t Ap t + Ap t Aq t ) 



(8) 



where Aq t := q t - n t (q) and Ap t := p t - Tr t (p). The 
covariance matrix V t changes in time deterministically 
according to the following Riccati differential equation: 

V t = AV t + V t A J + D 

- i(V t F T + hE T lm(C) r )(FV t + hlm{C)E), 



v a = v(p \y ), 



(9) 



where D := fi,ERe(C^C)£ T . Consequently, the optimal 
filter for the linear quantum system ^ is described by 
the closed set of equations Q and®, which is termed 
the quantum Kalman filter [Ilf/A remarkable fact 

is that the behavior of V t is determined without respect 
to the output Y t . This indicates that we can evaluate 
the quantum conditional expectation Vt = P( Pt \ yt ) by 
simply calculating the expectation Vt = (Pt)- Actually, 
as V t evolves deterministically, we can see 

(P t ) = (p{P t \y t )) = (v t ) =v t . 

Now, the quantum version of Linear Quadratic Gaus- 
sian (LQG) control problem is addressed as follows. For 
the linear quantum system driven by the quantum Gaus- 
sian noise, we aim to find an optimal control input it° pt , 
which is a function of the observations Y s (0 < s < t), 
such that the following quadratic cost function is mini- 
mized: 



J[ut] = ( J (^xjMx t + T -4)dt + ^x~TNx 7 



(10) 



The positive semidefinite matrices M > 0, N > and 
the scalar number r > reflect our control strategy; For 
example, if we are strongly restricted in the magnitude 
of the control input, a large value of r should be chosen. 
This problem can be solved by using the dynamic pro- 
gramming method. The optimal input is then given by 
u° pt = — (2/r)B T K t nt(x) , where the real symmetric ma- 
trix K t is a solution to the following Riccati differential 
equation: 

K t + K t A + A T K t - -K t BB T K t + \m = O, 
r 2 

K T = N. (11) 



Thus, we observe that the optimal control input u° pt is 
not a function of the entire observation history up to 
time t, but only depends on the solution to the Kalman 
filter ([7]) and (J9j) at time t. A controller that satisfies 
this desirable property is termed a separated controller. 
A general discussion on the optimality of the separated 
control is found in [271 ] . 



B. Quantum risk-sensitive observer 

The risk-sensitive control problem was originally for- 
mulated by Jacobson within the framework of the classi- 
cal control theory [l8| , and recently its quantum version 
was developed by James [25|, [2(| . The purpose is to de- 
sign an optimal control input such that the following cost 
function is minimized: 



J> t ] = (R T j T (e^)R 



Th 



(12) 



where R t is the solution to the operator differential 
equation dR t /dt = (fi/2)j t (a(ut))Rt and the parameter 
(i > represents the risk-sensitivity. The nonnegative 
self-adjoint system operators a(ut) and (3 are termed the 
running and terminal cost operators, respectively. In the 
classical case where the cost operators are scalar values, 
i.e., a(u t ) — a(u t ) and $ = f3, the cost function (fT2|) is 
reduced to 



JV[ Ut ) = ^exp^ J 



a(ut)dt + /i/3 



For this reason, Eq. \12\ is considered as a natural non- 
commutative generalization of the exponcntial-of- integral 
cost function. James has proved that the quantity (fT2"]) 
is expressed as 



cxp 



W 4(a)dt)n T {e^) 



(13) 



where E M denotes the expectation with respect to a cer- 
tain classical probability distribution (see [25|) and 7r^(») 
is a risk-dependent estimate of the system observable. 
The estimator is determined by the following equation: 

hdiTt(X) = ir?{£X)dt 

+ | [n?(Xa + &X) - 2^(X)^(a)] dt 

+ [^{Xc+^X)-4{X)4{c + ^)] 
x [dY t - K^c+c^dt]. (14) 

This differs from the filtering equation (j4]) in that the 
risk-dependent term is added to it. Therefore, tt^(X) is 
no longer the optimal estimate of the system observable. 
However, the risk-dependent term is indeed necessary in 
order for the cost function (|T2|) to be expressed only in 
terms of quantities defined on the system space that is 
driven by the output Y t . This implies that our knowledge 
about the system is tempered by purpose. 
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We now apply the above mentioned risk-sensitive con- 
trol theory to the linear system $Q with the following 
cost operators: 



1 



1 



a(ut) = —x T Mx H — u?, = —x T Nx. 
\ t) 2 2*2 

where M and N are 2x2 positive semidefinite matrices 
and r > 0. Then, through a lengthy calculation we obtain 
the corresponding observer equation as follows: 

di$(x) = (A + nVfM)i$(x)dt + Butdt 

VfF T + £ T Im(C) T ] [dY t - F^{x)dt), (15) 

where the time evolution of the symmetrized covariance 
matrix V/ 1 is given by 

vr = 



1 



V/ 1 A T + D 



(VfF T + hE T Im(C') T )(FVl l + film(C)S), 



a(W- t e t me), 



F " = P(P |%)- 



(16) 



Consequently, the risk-sensitive observer (JT4jl in linear 
case reduces to the closed set of equations, (IT51) and (ITBl . 
We also see that the cost function (TT51) is calculated as 



III. ROBUST OBSERVER FOR UNCERTAIN 
LINEAR QUANTUM SYSTEMS 

This paper deals with a linear quantum system such 
that specific uncertainties are included in the system 
Hamiltonian H and the system operator c as follows: 

H = \x T {G + AG t )x - x T Y.Bu u 
c= (C + AC t )x, 

where the real symmetric matrix AG* and the complex 
row vector AC t represent time- varying parametric uncer- 
tainties that satisfy the following bounds: 

(AG t f < gl, (18) 
(ReAC* t ) T (ReAG t ) < rj, (ImAG t ) T (ImAG t ) < r 2 I. 

(19) 

Here, the nonnegative scalar constants r%,T2, and g are 
known (/ denotes the 2x2 identity matrix) . By defining 



A A t := SAG t + Sim 



&AC t + ACjC + AC} AC t 



the dynamics of the system observable x t = [<jt Pt] T 
[jt(q) jt{p)Y ls represented as 



J >1 {ut}=W 



expl /i 



(^(x) T Mn?(x) + V -u 2 t )dt 



xexpU/ -Tr(M^)dt)4(e^) 



Note that the second integral in the above equation is 
a constant term as is deterministic. This is com- 
pletely a classical controller design problem and was al- 
ready solved by Jacobson [18j; the optimal control input 
that minimizes J^lut] is given by 



opt _ 



-B T KW(x), 



where satisfies the following Riccati differential equa- 
tion 

K? + K?A + A T Kt - -K?BB T K? + 



■2\iK\ 



tVfF 1 
h 



S T Im(C*) 1 



-FV^ + Im(C)S 



niKfVfM + MVfK?) = O, 
1 



K£ = - [(I - uNV^y-N + N(I - [iV^N)- 1 ] . (17) 

Therefore, it° pt is a separated controller composed of the 
solutions to the observer equation (fT5")) and the two cou- 
pled Riccati equations (fT6|) and (fTT)) . It is notable that 
these set of equations are identical to those in the quan- 
tum LQG optimal control problem when the risk param- 
eter /i is zero. In this sense, the LQG optimal controller is 
sometimes referred to as the linear risk-neutral controller. 



dx t = (A + AAt)x t dt + Bu t dt 

+ iE[{C + AC t ) T dB\ — (C + ACtYdBt] . (20) 

Moreover, the uncertainty is also included in the output 
equation as follows: 

dY t = (F + AF t )x t dt + dB t + dBj, 

AF t := AC t + AC*. (21) 

Here, we should remark that the drift and diffusion terms 
in Eq. ([20|) and the output equation (f2Tj) are affected by 
the common uncertainty ACt. This is because the quan- 
tum evolution is restricted to satisfy unitarity and the 
system matrices are thus strongly connected with each 
other. This is indeed an intrinsic feature of quantum 
systems that is not seen in general classical systems. 

Motivated from the structure of the Kalman filter 0, 
we aim to design a linear observer of the form 



dx t — Rxtdt + Bu t dt + kdY t . 



(22) 



where R and k are a matrix and a vector to be deter- 
mined such that the variance of the estimation error is 
guaranteed to be within a certain bound. The vector 
x t = [qt Pt] T G K 2 represents the estimate of the 
system observable Xt- Note that, as in the case of the 
risk-sensitive observer, Xt is not necessarily the optimal 
estimate of it- Furthermore, we here assume that the 
control input Ut is fixed to a linear function of the 
observer state, Ut — Lxt, where L is a row vector with 
the size 2. Then, an explicit form of (R, k) that enjoys 
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a guaranteed estimation error bound is provided in the 
following theorem. We remark again that the theorem 
can be easily generalized to the multi-particle case. 

Theorem 1. Suppose there exist positive scalars 
5i (i = 1,2) and e$ (i — 1, ... ,8) such that the following 
two coupled Riccati equations have positive definite 
solutions Pi > and P 2 > 0: 



(A + BL)P l +Pi{A + BL) T + PxQxPi 

+ D' + S 1 I = 0, 

A'P 2 + P 2 A' T + D' + S 2 I 
1 



(23) 



{P2F' T + MiS T Im(C , ) T )(F'P 2 + Milna(C)S) 

- P 2 (L T B T P 1 - 1 + P^BL)P 2 = O, (24) 

where the matrices A' and D 1 and the vector F 1 are de- 
fined by 

A' := A + (D + Q 2 + Q 3 )Pf \ D' :=D + Q 2 + Q 3 , 
F' := F + ^iIm(C , )SP 1 _1 . 

The definition of the matrices Qi (i = 1,2,3) are given 
in Appendix A: Eqs. ££H), jMJ, and ([A"6]l . The scalars 
/ii and \i 2 are given by \i\ — h + Ar\/e 2 and [i 2 = h + 
8ri/e 2 + he^, respectively. Then, the observer 



dx 



, = (A 1 - P 2 L T B T Pf x )x t df + Bu t dt 
1 r 



/'■2 



P 2 P' T + AiiS T Im(C , ) T (dY t -F'xtdt) (25) 



generates the estimate xt — [qt Pt] T that satisfies 

qt) 2 + (Pt - Pt) 2 ) <TrP 2 , (26) 
for all admissible uncertainties. 



lim ( (q t 

t—>oo 



Proof. We consider the augmented variable 

z t = [it it — Xt] T , where Xt and xt satisfy Eqs. 
([2U)l and (|2^|) . respectively. Then, z t obeys the following 
linear QSDE: 



dz t = {A+ AA t )z t dt + b A dBj 



b* A dB t 



(27) 



where 
A = 

b A = 



A + BL -BL 

A-R-kF R 

iS(C + AC t ) T 
iS((5 + AC t ) T - k 



,AA t 



AA t O 
AA t - kAF t O 



Let us now consider the symmetrized covariance matrix 
of z; V nm = (z„z m + z m z n }/2, (n,m = 1,...,4). This 
satisfies the following generalized uncertainty relation: 



[z t z t 



in. 



-E>0, E 



E E 
E E 



Noting B t $ = and the quantum Ito rule dB t dBj = hdt, 
the time evolution of V t is calculated as 



d_ 

7ft 



^-V t = (A + AA t )(z t zJ) + (z t zJ)(A + AA t ) T + hb A b T A 



(A + AAt) 

f> ih 
Vt + — 



(A + AAt) 1 



hb* A b A 



(A + AA t )V t + V t (A + AA t ) T + D + AD t . 



The matrices D and AD t are given by 
D 



' D 


D ' 




D 


D 


-h 



o 

km T km T 



mk T 

- m 



k T 



AD t = 



AD t AD t 
AD t AD t 



O 



kAml kAmJ 



kk 1 j 

Amtk T 

Am t k T 



where 

AD t ■= ftSRe 



&AC t + ACjC + AC J AC t 



>C}C + ACji 
m := E T Im(C') T , Am t := E T Im(AC' t ) T . 
Our goal is to design R and k such that the condition 
3X > 0, 

s.t. {A + AA t )X + X(A + AA t ) T + D + AD t < (28) 

is satisfied for all admissible uncertainties; in this case, 
it follows from the lemma shown in Appendix B that the 
relation limt^oo Vt < X is satisfied. For this purpose, we 
utilize the following matrix inequalities: For all X and 
the uncertain matrices satisfying Eqs. (fl8|) and (fl9|) . we 
have 



AA t X + XAAj < XQ X X 
AD t < Q 3 - 



(29) 
(30) 



The proof of the above inequalities and the definition of 
the matrices Qi (i = 1,2,3) are given in Appendix A. 
Therefore, the condition (|2"5]) holds for all admissible un- 
certainties if there exists a positive definite matrix X > 
such that the following Riccati inequality holds: 



AX + XA + XQ X X + D + Q 2 + Q 3 <0. 



Especially we here aim to find a solution of the form 
X = diag{Pi, P2} with Pi and P 2 denoting 2x2 positive 
definite matrices. Then, partitioning the 4x4 matrix ^ 
into ^> = (^ij) with 2x2 matrices \?.y, we obtain 

*n = (A + BL)P 1 + Pi (A + BL) T + PiQxPi + D', 
*2i = (A — R— kF)P 1 + D' - mkm J - P 2 L J B T , 
*22 = RP2 + P 2 P T + D' + /n(fcm T + mk T ) + ^i 2 kk T . 

Let us now assume that the Riccati equation ([23]) , which 
is equal to = —8x1 < 0, has a solution Pi > 0. Then, 
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the equality * 2 i = yields R = A' ~kF' - P 2 L T B T Pf 1 . 
Moreover, ^ 22 is then calculated as 



*22 = A'P 2 + P 2 A 
1 



/T 



M2 



/'2 



-D' 
Mi 

— m 



-PiF 



n 



1 



M2 

/i 1 m)(P 2 i ;l ' T 4 
- P 2 (L T B T P^ 1 + P^BL)P 



- — (P 2 P' T 
M2 



/'2 



Mi 
— '/ 

M2 



Hence, the optimal k that minimizes the maximum eigen- 
value of ^22 is given by 



k= — (P 2 F 

M2 



/T 



/11m = 



1 

/'2 



P 2 P' T + MiS T Im(C) T 



Then, the existence of a solution P2 > in Eq. 
directly implies ^ 22 = —S 2 I < 0. As a result, we obtain 
# = diag{— 5il, — S 2 I} < 0, which leads to the objective 
condition (|28|) . Therefore, according to the lemma in 
Appendix B, we have Hindoo V t < X . Then, as the 
third and fourth diagonal elements of the matrix Vt 
are respectively given by V33 = = ((q t — qt) 2 ) and 
V44 = (z\ 



((fit -Pt) ), we obtain Eq. (J26 



The basic idea to determine the form of the quantum 
robust observer (|25[) is found in several pap ers that deal 
with uncertain linear classical systems [2l[ I22I I23I [2^ |. 
However, the structure of the quantum robust observer 
differs substantially from that of the classical robust 
observer derived in [2l], [22|, [23|, [24[ ■ The reason for this is 
as follows. First, unlike the classical case, the covariance 
matrix Vt of the augmented system (|27|) , which is used to 
express the performance of the robust observer, must be 
symmetrized in order for Vt to be a physical observable. 
Second, the uncertainty AC* appears both in the drift 
matrix AAt and the diffusion matrix ADt in complicated 
ways; this is because, as has been previously mentioned, 
the system matrices are strongly connected with each 
other due to the unitarity of quantum evolution. The 
classical correspondence to the uncertain quantum 
system (f2TJl) and (j!?Tj) has not been studied. For this 
reason, the resulting robust observer ([25]) and the proof 
to derive it do not have classical analogues. Actually, for 
standard classical systems whose system matrices can 
be specified independently of one another, the process 
shown in Appendix A is unnecessary. 

We now present an important property that the 
quantum robust observer should satisfy: When the 
uncertainties are small or zero, the robust observer 
should be close or identical to the optimal quantum 
Kalman filter, respectively. This natural property is 
proved as follows. 

Proposition 2. Consider the case where the uncertain- 
ties converge to zero: AGt — > and AC t — ► 0. Then, 
there exist parameters Si (i = 1,2) and Ci (i = 1, . . . , 8) 
such that the robust observer ([25]) converges to the 
stationary Kalman filter ([7]) with Vt satisfying the 



Riccati equation V t — in Eq. ©. 

Proof. Let us consider the positive parameters 
€i (i = 1, . . . , 8) as follows: 

£ i = Vg-, £ 2 = max{^/rT, Vr-J}, £3 = max{^FT, x/ril 
e 4 = ri, e b = r 2 , e e = y/n, €7 = ^/7^, e 8 = yjr^. 

In this case, for example, the matrix Q\ is calculated as 

Ql = (Vff + max { \A7 , V^a} + r l + r-z) 1 
+ max{VFT, V^KCjCi + CjC 2 ), 

which becomes zero as g — > 0, r± — > 0, and r 2 — * 0. 
Similarly, in these limits, we have Q 2 — ► 0, Q3 — * 0, jix — > 
H, and /12 — * h. Then, since Eq. (|23|) is equivalently 
written as 

Pi\A + BL) + {A + PL) T Pf 1 + Q a 
+ P 1 -\D' + ( 5 1 /)Pf 1 -O, 

the limit Qi ^0 implies that the solution of the 
above equation satisfies Pf — > 0. We then obtain 
A' —fA,F'—f F, and D' -> D. Therefore, in this 
case, Eq. (IM1) with 52 — is identical to the Riccati 
equation Vt = in Eq. @. The robust observer (f2l)|) 
then converges to the stationary Kalman filter (JT)) with 

v t = p 2 . m 

The above proposition also states that we can 
find the parameters 8i and such that the robust 
observer l|25p approximates the stationary Kalman filter 
when the uncertainties are small, because the solutions 
of the Riccati equations ([23]) and ([24]) are continuous 
with respect to the above parameters. 

We lastly remark on the controller design. In Theo- 
rem 1, we have assumed that the control input is a linear 
function ut = Lx t . This is a reasonable assumption in 
view of the case of the LQG and risk-sensitive optimal 
controllers. Hence, it is significant to study the optimiza- 
tion problems of the vector L such that some additional 
specifications are further achieved. For example, L opt 
that minimizes the upper bound of the estimation error, 
Tr P2, is highly desirable. However, it is difficult to solve 
this problem, since the observer dynamics depends on L 
in a rather complicated manner. Therefore, the solution 
to this problem is beyond the scope of this paper. 



IV. EXAMPLE— FEEDBACK COOLING OF 
PARTICLE MOTION 

The main purpose of this section is to show that there 
actually exists an uncertain quantum system such that 
both the robust observer and the risk-sensitive observer 
perform more effectively than the Kalman filter, which is 
no longer optimum for uncertain systems. Moreover, we 
will carry out a detailed comparison of the above three 
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observers by considering each estimation error. This is 
certainly significant from a practical viewpoint. 

First, let us describe the system. The control objective 
is to stabilize the particle position q at the origin by con- 
tinuous monitoring and control. In other words, we aim 
to achieve (7r t (g)) = (q) — with a small error variance. 
The system observable is thus given by 

c = q, i.e., £=[10]. 

For the Hamiltonian part, H = H llcc + H contro \ we as- 
sume the following: The control Hamiltonian is propor- 
tional to the position operator: 



H 



control 



-u t q, 



i.e. 



B = 



(31) 



where u t — Lxt is the input, and the free Hamiltonian is 
of the form H hco — 2p 2 + V(q), where V(q) denotes the 
potential energy of the particle. In general, the potential 
energy can assume a complicated structure. For example, 
Doherty et al. [28| have considered a nonlinear feedback 
control problem of a particle in a double-well potential 
V(q) = q* — q 2 . Since the present paper deals with only 
linear quantum systems, we approximate V(q) to the sec- 
ond order around the origin and consider a spatially local 
control of the particle. In particular, we examine the fol- 
lowing two approximated free Hamiltonians: 



free 



2p 2 -0.05<f, H. 



free 



2p 2 + 0.05<f 



The former is sometimes referred to as an anti-harmonic 
oscillator, while the latter is a standard harmonic oscilla- 
tor approximation. The system matrices corresponding 
to H frcc are respectively given by 



Gi 



-0.05 
2 



Go 



0.05 
2 



In the case of the harmonic oscillator Hamiltonian, the 
system is autonomously stable at the origin. In contrast, 
in the case of the anti-harmonic oscillator, the system 
becomes unstable when we do not invoke any control. 
However, it is observed that the control Hamiltonian ([31]) 
with an appropriate control input can stabilize the sys- 
tem. An example is the LQG optimal controller with the 
following tuning parameters of the cost function CH 



M = 



3 
1 



1 

r=-, N = 
5 



2 




(32) 



Figure I illustrates an estimate of the particle position in 
both the unstable autonomous trajectory and the con- 
trolled stable trajectory; in the latter case, the control 
objective (^t(q)) — is actually satisfied. 

Second, we describe the uncertainty included in the 
system. In particular, we consider two situations in which 
uncertain Hamiltonians AHi = — \fd~ t q 2 and AH2 = 
\/dtq 2 are added to Hi and H2, respectively. The un- 
known time-varying parameter dt is bounded by the 




FIG. 1: An example of the unstable autonomous trajectory 
(dot line) and the controlled stable trajectory (solid line) 
shown by n t (q). 



known constant g > 0, i.e., d t G [0,g]. Regarding 
the uncertainty in the system operator c, on the other 
hand, we assume ACt — 0, Vt. In this case, we can set 
Qi = (\I,Q2 = {g/e\)I, and Q3 = by choosing the 
parameters shown in the proof of Proposition 2. 

The comparison of the three observers is performed 
based on the following evaluation. For the Kalman filter 
and the risk-sensitive observer, we evaluate the station- 
ary mean square error between the "true" system and 
the estimator for the "nominal" system corresponding 
to dt = (see Appendix C). In both cases, the tun- 
ing parameters in the cost function are set to Eq. (|32p . 
Next, for the robust observer, we evaluate the guaran- 
teed upper bound of the estimation error Tr P2 in Eq. 
([26]) . The control input in the robust observer is set to 
the stationary LQG controller for the nominal system: 
Ut = Lx t = — (2/r)B T KooXt, where Koo is the station- 
ary solution of Eq. (fTTj) . 

Let us now describe the simulation results. First, we 
consider the case in which the total system Hamiltonian 
is given by H = fff oc + H contro1 + AH X . Table I lists 
the three estimation errors mentioned above for several 
values of g. Here, the uncertainty d t is set to the "worst 
case" dt = g for each value of g. In the first row of the ta- 
ble, the notation "N/A" indicates that the solution of the 
Lyapunov equation (|C1[) does not satisfy W + iKS/2 > 0. 
This implies that the error dynamics between the un- 
certain actual system and the nominal Kalman filter is 
unstable. In other words, the Kalman filter fails in the 
estimation. It should be noted that two excessively large 
values of the estimation error, which appear in the first 
and second rows, indicate that the error dynamics is 
nearly unstable. Therefore, it can be concluded that 
the Kalman filter and the risk-sensitive observer for the 
nominal system do not work well when the uncertainty 
dt (= g) assumes a large. On the other hand, as shown in 
the third row in Table I, the robust observer is not very 
sensitive to the magnitude of the uncertainty and pro- 
vides a good estimation even when g is large. The above 
discussion suggests that the robust observer is possibly 
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TABLE I: Comparison of the Kalman, risk-sensitive, and ro- 
bust observers, denoted by "KAL", "RSK", and "ROB", re- 
spectively. The free Hamiltonian of the system is approxi- 
mated by the anti-harmonic oscillator. In order to calculate 
the guaranteed upper bound of the estimation error of the 
robust observer, Tr P%, parameters 5i and 82 are fixed to 0.1, 
and ei is selected such that Tr P2 takes the minimum value. 
The risk-sensitive parameter is fi — 0.3, and the Planck con- 
stant is set to unity: h — 1. Note that both the robust ob- 
server and the risk-sensitive observer are not identical to the 
Kalman filter even when g — 0, because the parameters 62 
and /1 are now set to non-zero values. 



g 0.00 


0.20 


0.38 


0.60 


0.80 


0.97 


KAL 1.43 


2.38 


40.88 


N/A 


N/A 


N/A 


RSK 1.48 


1.82 


2.21 


3.19 


6.07 


61.27 


ROB 1.73 


3.32 


4.74 


7.04 


10.12 


14.13 


TABLE II: Comparison of the three types of estimators 


in the 


case of the harmonic oscillator Hamiltonian 


. All parameters 


of the estimators 


are set to the same values 


in Table I. 




g 0.00 


0.20 


0.40 


0.60 


0.80 


1.00 


KAL 1.40 


1.37 


1.40 


1.44 


1.47 


1.50 


RSK 1.44 


1.38 


1.38 


1.39 


1.40 


1.41 


ROB 1.68 


3.23 


4.79 


6.84 


9.80 


14.48 



the best option for dealing with a large uncertainty. In 
other cases, the risk-sensitive observer should be used. 

Next, we consider the second example, in which the 
total system Hamiltonian is given by the harmonic os- 
cillator H = H$ ee + H coat101 + AH 2 . In this case, it is 
immediately observed in Table II that the estimation er- 
rors of the robust observer are always greater than those 
of the others, while the risk-sensitive observer shows a 
good performance, particularly when g assumes a large 
value. Hence, in this case the risk-sensitive observer is 
the most appropriate. 

An interesting feature of the robust observer is that in 
the case of both the harmonic and anti-harmonic Hamil- 
tonians, it provides almost the same trend in the estima- 
tion errors with respect to g, whereas the Kalman filter 
and the risk-sensitive observer produce drastically differ- 
ent trends in the errors. This indicates that the structure 
of the robust observer is designed such that the estima- 
tion error is insensitive to the stability property of the 
system. However, this design policy sometimes leads to 
the over conservative stability of the error dynamics, and 
the estimation performance eventually reduces. 

V. CONCLUSION 

In this paper, we have considered a linear quantum 
system subjected to time-varying norm-bounded para- 
metric uncertainties and developed a quantum version of 
the robust observer. Although in the linear case much 
of classical control theory can be applied to quantum 



systems, due to the unitarity of quantum evolution, the 
quantum uncertain system must have a specific structure 
with respect to the uncertainties, and its classical corre- 
spondence has not been studied; the resulting quantum 
robust observer has thus no classical analogue. The ob- 
server differs from both the optimal Kalman filter and the 
risk-sensitive observer; however, it guarantees the upper 
bound of the variance of the estimation error. We then 
investigated the robustness property of the three esti- 
mators mentioned above by considering a typical quan- 
tum control problem — feedback cooling of particle mo- 
tion. This examination clarified that the robust observer 
is superior to the others when the autonomous system 
is unstable and is subjected to an unknown perturba- 
tion with a large magnitude. Therefore, we can conclude 
that the robust filtering method originally developed for 
classical systems is actually very effective for quantum 
systems as well. This fact implies that several robust 
control techniques in classical control theory (e.g., (29J) 
will be applicable to uncertain quantum systems. 
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APPENDIX A: PROOF OF EQS. ((29]) AND ([30)) 

At first, we derive a simple yet useful matrix inequality. 
For any real matrices X and Y, we obviously have 

(V-eX-±Y) T (V-eX-±Y)>0, 

where e > is a free parameter. The above inequality 
immediately leads to 

X T Y + Y T X < eX T X + -Y T Y. (Al) 

e 

Next, let us define 

C*i = ReC, C 2 = ImC*, AC\ = ReAC, AC* 2 = ImAC*. 

(In this appendix, we omit the suffix t for simplicity.) 
Then, the conditions (fl"9|) are represented by AC J AC i < 
ril (i = 1, 2). Note that they lead to the scalar inequali- 
ties: ACiACj <n {i = 1,2). 

Now we are at the point to prove. Let us first derive 
the inequality ([2^]) . By a straightforward calculation, we 
obtain 



AA = T,AGE + OiAJiS + SA J 2 6 2 

+ tACjAC 2 E - EACjACiS, 



9 



where E = -[E E] e R 4x2 , £=[/0]e 



p2x4 



and 



The matrices Qi and Qi are defined as follows. 



0i = 
AJi = 



EC?" -EC 
A<7 2 



2fc 



e 2 = 



c 2 o T 

-Ci T 



AJ 2 = [ACj ACj]. 



We here denoted T = [0 0]. Accordingly, the matrix 
AAX + XAA T is now represented by 

AAX + XAA T 

= (tAGE)X + X{tAGE) T 

+ (GiAJiE)Z + X(Q 1 AJ 1 E) T 

+ (tAJ 2 Q 2 )X + A > (EAJ 2 9 2 ) T 

+ {tACjAC 2 E)X + X{tACjAC 2 E) T 

- (EACj ACiE)X - X{tAC 2 T AC\E) T . (A2) 

We are then able to apply Eq. (|A1[) to evaluate bounds of 
each line in the above equation. For example, the second 
line has the following bound: 

(T,AGE)X + X(T,AGE) T < ^XE 1 EX + — EAG 2 E T 

£l 

< exXE^EX + — EE T , 
£i 

where here the assumption on the uncertainty (|18j) was 
used. The free parameter t\ > should be tuned appro- 
priately. Next, for evaluating the third line of Eq. (IA2[) . 
we remark the following: 



AJxAjJ < diag{2r 2 , 2n}. 



(A3) 



This inequality is easily seen; the relations AQACj < 
Ti (i = 1, 2) lead to 

det diag{2r 2 , 2n} - A.hA.jJ 

= ||ACi|| a ||ACfe|| a - (ACi, AC 2 ) > 0. 

Here, ||»|| 2 and (•, •) denote the standard Euclidean norm 
and inner product, respectively. By using Eqs. (|A1|) and 
13]), we then obtain the following inequality: 

(6iAJiE)J? + ^(eiAJiE) 1 " 



< e 2 XE 1 EX + — 9iA Ji A J^e]" 

£2 

< e 2 XE T EX + — 9idiag{r 2 , n}Qj. 

£2 



For the other lines of Eq. (IA2[) . we can use the same 
manner to have their bounds that do not depend on the 
uncertainties. As a result, we obtain the objective in- 
equality A AX + XAA T < XQ 1 X + Q 2 , where 



Qx ■= 

Q 2 := 



Qx O 
O O 

Q2 Q 2 
Q2 Q 2 



An 

£2 



O 

km T kr 



mk T 



2/c/fc 1 



Qx := (ex + e 2 + e 4 + £5)/ + e 3 (CjC 1 + CjC 2 ), (A4) 



Q2 := 



?"l + ?"2 ^1^2 



£3 



£4 



— T,(r 2 CjCi +riC 2 C 2 )T, T . (A5) 
£2 



The parameters e, > (i — 1, . . . , 5) should be chosen 
appropriately. 

Let us next derive Eq. (|30p . Similar to the previous 
case, we use Eq. (IA1[) to obtain a bound that does not 
depend on the uncertainty. First, we immediately obtain 



AD = 



< /iE 



C^ACi + ACjC! + CjAC 2 + ACjC 2 

+ ACjAC! + ACjAC 2 ~\ E T 

eeCjCx + + e 7 CjC 2 + —I 
£6 £7 



+ (n+r 2 )I E T =: Q' 3 



!3- 



where cq > and €7 > are free parameters. This readily 
leads to 



AD AD 
AD AD 



< 



Q'3 Q'3 
Q'3 Q'3 



Also, setting X = [0 T k T ] and Y = [-Am T - Arn T ] in 
Eq. (|Alj) . we obtain the following inequality: 

Amk T 
kAm T kAm T + Amk T 

AmAmJ /e$ AmAm T /es 
AmAmJ/eg AmAm T /eg + e%kk T 

r 2 I/e s r 2 I/e$ 
r 2 I/e s r 2 I/e s + e s kk T 



< h 

< h 



where es > 0. Consequently, AD is bounded by AD < 
Q3, where 



Q3 ■= 



^>3 W3 
Q3 Q3 



he* 



O O 
O kk T 



and 



Q3 :=n(^ + ^ + ^+rx+r 2 )l 

+ KEieeCjC! + e 7 (7 2 T C' 2 )E T . 

APPENDIX B: UPPER BOUND LEMMA 



(A6) 



We consider a matrix-valued differential equation of 
the form 



P t = AP t +P t A T + BB T . 



(Bl) 



Lemma. Suppose there exists a positive definite matrix 
X > such that the inequality AX + XA T + BB T < 
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holds. Then, Eq. (|Bip has a unique stationary solution 
that satisfies 

lim P t < X. 

t — >oo 

Proof. We readily see that the matrix A is strictly sta- 
ble; any eigenvalue of A has a negative real part. Now, let 
us define SP t := X — P t . Then, by using the assumption 
we have 



SP t = -(AX + XA J + BB T ) 
> A8P t + 8P t A T , 



^A't 



A5P t + SP t A T 



We then obtain 



which yields 5Pt > e SPoe 
limt^oo SPt > since A is strictly stable. This shows 
the assertion. ■ 



APPENDIX C: NOMINAL-TRUE SYSTEMS 
DIFFERENCE 

The objective here is to characterize the stationary 
mean square error between the "true" system ([20)1 and 
the risk-sensitive observer (|15p specifically designed for 
the "nominal" system ([5]). For this purpose, we calculate 
the symmetrized covariance matrix of the error vector 
e t := x t — 7rf (£), where Xt and 7rf (x) are generated from 
Eqs. ((20)) and (fT5|) . respectively. Particularly, we now 
focus on the stationary observer. Thus, let us assume 
that the two Riccati equations p^|) and (TlT)) have unique 
steady solutions and K^, respectively. Then, defin- 
ing 

K := V^F T + £ T (ImC<) T , L Q := ——B T K£ D , 
h r 



the stationary risk-sensitive observer is described by 

dic?(x) = {A + nV£M + BL )^{x)dt 
+ b (dY t -F^(x)dt). 

We then see that the augmented vector £ t = [xt et) 
satisfies d( t = A (tdt + b dB t + b dBj, where 

A n 



bo = 



A + AA + BL -BL 
AA-nV£M A + f iV£M-b F 



iSC 1 



Let Vt be the symmetrized covariance matrix of £ t . As 
mentioned in the proof of Theorem 1, this matrix satisfies 
(CtCt ) = Vt + ihfj/2. By using this relation, we obtain 
dVt/dt = A V t + V t A T + D , where D a is given by 



D„ 



D D~ 






D D 






O 




[Mm(C)£] T 


6 Im(C)S 


Mm(C)£ + [6 Im(C)S] T - b b T \ 



As a result, the variance of the estimation error is given 
by lim t ^ OQ (e7e t ) = W33+W44, where W is the stationary 
solution of the following Lyapunov equation: 



A D W + WAl +D = 0. 



(CI) 



The estimation error between the true system and the 
Kalman filter designed for the nominal system is immedi- 
ately evaluated by setting /1 = in the above discussion. 
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